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: 3 The fields of Missile Ballistics and Missile Dynamics are so 
closelly relisted to each other that one has to define both fields 
carefully, if one wants to draw a dividing line between them. 

oe Both disciplines deal with the motion of missiles in space. (The 

ue. word space is used here in the sense customary in mathematics, not 

” qs space t“avelsrs use it, This means our space may be filled by 
eC air of anyfdensity,, including zero density.) 


ee) 
ee Dynamics, meaning motion caused by acting forces, could include 
afiz, m ie of missiles. Ballistics, derived from the Greek word 
Aol) Ip t, to throw, could also be applied to any phenomenon 
94 with the hurling of a missile through space. Sticking to 
Pissed meaning of inese words we define: 


fissile Ballistics is the theory of the trajectories of 
pat foaghy th of ais: eeee: ~ Data of missile ballistics 


ge = Missile mamic3 is the iiaase of the fine-structure of 
Bei trajectories dnoludine rotary motion of the missile about, 
“itg er ofvgravity and elastic deformations. -As a rule, data 
2 dynamics can be observed only by instruments installed 
sthe missile. 


mrity regulations have separated the science of ballistics 
om other sciernses a long time before airplanes and missiles even 


h ; fro 

hs > jexisted. Therefore ballisticians had to operate their own test 
aaa bs facilities and did not see any need for adapting their terminology 
FES and symbols to neighbouring fields, aerodynamics, for instance, when 
CM 7 these fields were buiit up. We will disregard this fact and will use 


the terminology of conventional aerodynamics, where they apply. 


=i ghed information will be given in this course, 


1cs in general is divided in three fields: 


Interior Ballistics 
Exterior Ballistics 
“4¢ End Ballistics 


Interior’ B.llistics deals with the phenomena which oceur inside the 
gur barrel, which accelerates the projectile. End Ballistics treats 
the effect of the projectile on its surroundings at the end of the 
trajectory. This course is restricted to the Exterior Ballistics 

of Missiles. 


Extericr ballistics of projectiles usually includes the 
aerodynamics of the shell, as for instance, drag problems and the 
selection of an optimum shape of the shell. In missile ballistics 
these problems are treated under the heading of missile aerodynamics. 


A remark has to be made on BALLISTIC TRAJECTORIES, which are 
defined as trajectories without aerodynamic sideforces, or since 
missiles with ballistic trajectories usually have two planes of 
symmetry, with zero angle of attack and zero angle of sideslip. 

Ou Missile Ballistics will not be limited to this kind of trajec- 
tories. If one looks closely into the forces acting on this type of 
trajectories, one would find that all types of trajectories which 
have curvature or are subject to natural wind, have lift forces. 
This is demonstrated by the 
sketch on the side. The angle 

of attack varies along the 
longitudinal axis of the missile, 
because the angular velocity, Gy » 
about the lateral axis causes a 
local change of the angle of 
attack. Consequently there is a 
pressure distribution along the 
missile, which has a resultant 
side-force, which is different 
from zero. 


aN 

ANCLE OF ATTACK 

Fig, 1. The angular 
velocity about the 
lateral axis produces 
a lift force. 


THE ROCKET EQUATION 


A rocket is a reaction type engine which produces thrust, T, by 
ejecting gases of burned propellant with an ejection speed, ce The 
rocket contains all ingredients necessary for combustion in its 
propellant and is in every respect independent from the presence of 
air. 
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We derive the equation of motion of a rocket in a space free 
from air and gravity from the principle that in this type of space 
the center of gravity of a system retains its location in space as 
long as internal forces act only. If one considers the expelled 
gases part of the system, the rocket thrust is an internal force. 


z 


Fig. 2. Snapshot of a rocket and its trail of expelled 
gases at the time, t. 


At the time, t, the rocket has traveled a distance, s. In 


each second the mass ‘ 
Mon =m [Bese] 
Y 


has been ejected with a speed, c, relative to the rocket. The gases 
which were ejected at time, t=r-, have a speed C— Aj lt) to the 

right relative to the raat, Therefore they have traveled at the 
time, t, a distance d= (c~- A le)) [t~ t) 


The condition, that the center of gravity of the system has 
not, moved since the beginning of the motion can be expressed by 
the statement: the moment of the masses about the center of gravity 
location at the beginning of the motion is zero. This is expressed 
in a farmula: 


boy pi . 
m b)edl4) = - [m de] (e-dery|/t-7) - ale) 
0 


Differentiation with respect to t gives, if one considers that m 
is a constant and that the variable, t, 1s oe contained in the 
upper limit of the integral: 


BA rma = gb - [i deen se 


A second differentiation with respect to eh ctetdas 
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vat Md = - Mo + WA 


(ayn = mh Se as She “ .. ~ = We ms we 
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If one writes for s the variable v, one obtains the rocket equation 


(1) n> ne ates 


This equation can easily be integrated by separation of the vari- 
ables 


or <= 


Usually the exhaust velocity, c, is not given, but the 
specific impulse, I. The name of the specific impulse may need 
an explanation. The specific impulse has not the dimension of 
an impulse, which would be [1b e sec} but rather the dimension of 
time, The physical meaning of specific impulse is the length of 
time one pound of propellant can produce one pound of thrust. 


The relation between the specific impulse and the exhaust 
velocity, c, can be obtained, if one expresses the thrust in two 
ways, by means of c and by means of I. 


The only force acting on the rocket under consideration is 
the thrust. According to Newton's law of motion is 


(mass) + (acceleration) = acting forces. 


Consequently, the term -cm in equation (1) must be the thrust, 
T. The impulse produced by the thrust in the time, At, is 
T° At. During this time the weight of the propellant has 
changed by -mg 4t, which contained the impulse -Img At. 
Both impulses are equal and one obtains by comparison of both 
expressions? 
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g, here, is the constant of gravitation, which is on the earth 
approximately 


(5) gn ato firseet 


When all propellant has been burned, the mass of the rocket 
has decreased to the value m,, the empty mass; By that time the 
velocity has increased to the end velocity, v,. Equation (2) 
gives the relation between the so called Massratio, m,/m, and 
the end velocity, v, : 


Ne 
If one instroduces the mass of the propellant, m,,; one has 


and instead of (6) 


(6a) UEP MATEO t edne (it P ) 


If one has a two-stage rocket, the end speed of the first 
stage, (ve)y5 becomes the initial velocity, (v,)o of the second 
stage. One obtains: 


mag t Oy eA [e) Ae) 7 


The product (my /mg), * (mg /me), is called the compound mass 
ratio of the two stage rocket. 


It is obvious, how equation (7) has to be generalized, if 
one wants to use it for a rocket of more than two stages and for 
configurations, which have different exhaust velocities in the 
different stages. 
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TRAJECTORIES WITHOUT ATR FORCES AND WITHOUT THRUST 


The equations of the trajectories contain a number of experi- 
mental functions and can therefore be solved numerically only. 
Before the development of the electronic computers, the ballistic 
offices employed a great number of people who had to operate desk 
computers in order to produce the necessary numerical solutions. 
It is obvious that in a course on missile ballistics of only 150 
minutes, there is no time for discussing the many different 
approaches for the numerical solution of the trajectory equations. 
On the other hand, it is the purpose of this course to give you a 
bird's eye view of the different methods used in missile ballistics 
and to provide you with the tools which enable you to solve some 
practical problems and to understand the more involved problems. 
For these reasons we start out here with the discussion of the flight 
of a ballistic missile in empty space after burnout, which covers the 
longest part of the trajectory of a ballistic long range missile. 
In the past two different methods have been developed and used for 
the calculation of this type of trajectory. 


(a) The original approach was the so called CORRECTION METHOD. 
The trajectory is first calculated under the assumption that the 
gravity force is constant with respect to magnitude and direction. 
The resulting trajectory is then a parabola. Afterwards corrections 
are applied which take care of the variation of the gravity force. 


(b) The ORBIT METHOD treats the missile like a satellite fly- 
ing in the gravitational field of the earth. 


It is obvious that the corrections on the parabola trajectory 
beGome the larger the longer the range of the missile. Experience 
has shown that the second method works faster than the first one as 
soon as the range of the missile exceeds some hundred miles. 
Therefore we will discuss the ORBIT METHOD only. 


We shall approach the discussion of this method in two steps: 
First it is assumed that the earth does not rotate, later the 
rotation of the earth is considered. 


ORBIT FORMULAE 


It is assumed that the part of the trajectory we deal with here 
begins in a point, B, which is located high enough above sea level 
so that one is allowed to neglect aerodynamic forces in this altitude. 
The velocity in point B may be given according to magnitude and 
direction. 
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Kepler's first law states that the trajectory is an ellipse 
(Ref. 9) with the center of the earth located in that focus of the 
ellipse, which is further away from the trajectory part of the 
ellipse. (See Figure 3) The rotation of the earth has no effect 
on the form of the trajectory because there is no connection with 
the earth after the missile has left the atmosphere. The only 
effect caused by the rotation of the earth is that one would measure 
from the rotating earth a different initial velocity. This correc- 
tion is given later in explicit forn. 


All symbols used in the following formulae are defined in Fig. 3 
and Fig. 4. It should be noted that the beginning, B, and the end, T, 
of the part of the trajectory, for which the following formulae are 
valid, have the same altitude above ground. 


The trajectory equation in polar coordinates is 


2 
(8) res hl 


g, (1- =cos ee 


8, 1s the gravitational acceleration in the initial point, B. It 
suould be noted that the gravitational acceleration decreases 
according to the formula 


| 2 
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where r, is the radius of the earth and go the gravitational accelera- 
tion at sea level. 


Introducing the abbreviation 


(10) J *ra/ "nl 


for the ratio of the components of the velocity at B and the geodetic 
parameter 


2 
(11) k= M4. T 40 
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the numerical eccentricity,E , of the trajectory ellipse can be 
written 


2 ate 
(12) Siri ld 2dedmekrd soc 
and the half angles, > of the trajectory arc of the ellipse 
(13) tg p = Arai 
1-k 


Equation (13) gives the range, R, of the trajectory 


The radial and the horizontal velocity, in a point with the centri 
angles, >» are respectively 


(15) ve = dr/dt “Vil = sin p 
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If one introduces in (8) the geodetic parameter, k, one has 


(16) Vn 


i 
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(17) rosfy k 
1-€ OEah 


and the peak altitude for ¥ = O becomes 


(18) Tmax “TL = Ty ( k -l) 


The speed in the peak point is 


(19) Mie Me ee: 


hp hl 


The flight time from the initial point, B, to the target point, 
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SOLUTION OF THE TWO BASIC BALLISTIC PROBLEMS 


We introduce spherical coordinates on the sphere, which is 
concentric with the earth and contains the initial point, B, and 
the target, T, (See Figure 4). 


DP = geographical latitude 


C 


geographical longitude, zero for the meridian of 
the initial point B. 


if 


V angle defining the direction of the flight path 
in the initial point. 
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The FIRST BASIC BALLISTIC PROBLEM an then be formulated as 
follows: 


Given are the coordinates r , 01; and Q, = O of the initial 
point, B, the velocity components v,; and v,, in point B, and the 
heading direction, VY; of the speed v] . Wanted are the coordinates 


D > ie OF +4D and Q, -4© of the target point, T. 


The solution for not rotating earth is given by the orbit 
formulae. One calculates the parameters ¢ and K by means of 
equations (10) and (11) and obtains the range angle 2% , from (13). 
The numberical eccentricity follows from (12) and the flight time, 
te , from (20). 


The coordinates of the target point, T, can be derived by simple 
spherical trigonometry from the triangle BIN in Figure 4. One has 
the relations 


(21) sin (2,440) - sin ?, cos 24, + cos P > sin 27; sin 


(22) sin A@ - Sin2 WY» cos 
E cos (P +a ) 


These are the formulae for not rotating earth, where the velocity 
components v;) and vj, ] measured from the earth are the same as in 
inertial space. In ie following the effect of the rotation of the 
earth is considered and the velocity components v,) and Vl? which are 
now measured from the rotating earth, are converted into the components 
Vr, and Vhg relative to the not rotating space. This conversion is 
obtained by considering the circumferential speed 


(23) uzr, Lb) cos y) 
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of the earth in the initial point, B. 4 is the angular velocity of 
the earth, whichhas the numerical value (= 0.72905 1074 rad/sec. 
One obtains 


(2h) bbe — igh 
Rage Trae Al A 
(25) = + 
Ving Sans LF2 Mcosy HM 
where 
Cg gece 
Hot Vhy 4 
The direction, , has also to be corrected for earth rotation 


and converted into D3 which may be determined by either of the 
following formulae 


me a0 sin y 
(27) sin Y= // + dpe cory + 


t+ COS 


(28) cos | 
Po indem cosy 47 


Likewise the flight parameters of and k are affected by earth 
rotation. They have to be cay into 


(29) 1, = /1 pad Mm Cosy tat 


| 3] 
(30) Ae ane Oa phe cosy te") ae 


Substituting (29) and (30) into (12), (13), and (20) for o and k, 
one obtains the eccentricity, €, , the range angle, 2Y, , and 
the flight time, tr , of the trajectory in the nonrotating space. 
(21) and (22) provide the coordinates on the sphere through B and T. 


While the missile moved from Bto T, the earth has rotated in © 
underneath the trajectory by an angle (): tro: Consequently the 


coordinates of the actual target point, T, on the rotating earth are 
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If one wants to know the errors caused by neglecting the 
rotation of the earth, one has to do it the hard way performing the 
calculations for rotating and for not rotating earth. There is no 
easier way: | 


FLIGHT TIME AND RANGE ERRORS OF PARABOLA TRAJECTORIES 


The errors caused by simplifying the orbit trajectory so that it 
becomes a parabola can be given in relatively simple form for the 
flight time and for the range. The flight time of a parabola trajectory 
is 2Vn1/g1° Comparing this value with (20) one sees that the flight 


time is in error by the factor 


Jj}—-s2 iL 
Lk on Bike tan 3) s ly gz 
J(1~s*) fiz l-E4 (ltée) tanh 


The range of a parabola trajectory is (2v 4 Vnp/81- Tr one 
deals with a missile for which the range angledy, is small enough so 
that tan y, may be replaced by Y one obtains’ from (13) thes 


pa a 
Ra 2 vary = “rl hl 
g (1-k) 


Consequently the range of the parabola is wrong by a factor 1/(1-k). 
THE SECOND BASIC PROBLEM OF BALLISTICS! 


The basic ballistic problem occurring in practice is an inversion 
of the first ballistic problem discussed so far. It is called the 
second ballistic problem and can be formulated as follows: 


The coordinates of the initial point of the trajectory and the 
coordinates of the target point are given. This means one knows 
rl; rae: A Q, A4@. Wanted are the velocity vector and its heading 
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in the initial point of the trajectory, this means v,), V,, and 
4 have to be determined. 


There is no closed solution known to this problem. One 
proceeds as follows. The rotation of the earth is first neglected 
and one calculates w and ¥, from (21) and (22) and from (13) a 
pair of corresponding values and k. The necessary velocity 
components are then given by (10) and (11). Now, (20) gives the 
flight time and one can correct the results for earth rotation. 
After this, the calculation is repeated. In practice, one calcu- 
lates large families of solutions of the FIRST BALLISTIC PROBLEM and 
selects from these solutions: the one which has the desired initial 
point and target point. 


OPTIMUM RANGE 


It is of interest to know, under which elevation angle one has 
to start the trajectory, if ore wants to obtain the longest possible 
range with a given initial speed. The answer can be derived from (13). 
Introducing the angle Vs between the vertical and the velocity vector 
and substituting in (13) 


“hy = Vo sin a 


V 
=) 


Vg cos wr, 


and differentiating with respect to ve one obtains for sition’ 77, 


the condition yy: 
/ Vo 
(32) serierady) Secligs F eenetras LUN Bee Pa ihe eres yee, eee NET 
Anes 9, ef 
x : Ldy 
are 
B) Oo 
For missiles, for which _‘o_. is. small the optimum elevation 


2rig 
angle of the vector of the initial speed, 90° - qn, » is only 
somewhat smaller than 45°. For missiles with longer range the 
optimum elevation angle can become appreciably smaller than 45°. 


NUMERICAL EXAMPLE 


The formulae may be illustrated by some numbers taken from Ref. '7. 
The example is a certain variation of the V-2 missile. 
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Initial point is the burnout point in 30 km altitude, which 
gives 


Es = 6 9 400 km 
components of burnouts velocity \y 
SAA iw’ y 
Sea ie 
Yr, = hy = 1,000 m/sec 
g, = 9.6963 m/sec 


D = 50° 


j ” 
For a parabolic trajectory is the range 


R= (2v,, Vi /€1 = 206.264 km 


The side deviation depends upon the flight direction, which was 
symbolized by WY . It lies between 3.64% and -1.3%, the range 
deviations between 2.81% and 0.24%. 


The flight time is for parabolic trajectory 


te - 2v.,/ g, = 206.264 sec (same number as range 
f yi oS) 
pied P because v,, = 1,000) 


Without rotation of the earth the flight time is for orbit 
trajectory “ait 
pepe tee = 211.935°see 


\ 


If the rotation of the earth is considered, the flight time depends 


upon the direction of flight and varies between 210.170 and 214.363 sec. 


CORRECTIONS FOR AIR DRAG 


For the ascent part in dense air(from the ground to point B in 


Figure 3) the influence of air drag and lift has to be taken care of in detarl 
as will be seen later. For the reentry part of the trajectory (from T 

to tae ground) the same is true. There exist, however, cases where the 
burnout point, B, is not quite high enough go that one would not like to 
neglect the air drag completel For such cases the following formulae 
are useful. (Fhoiegh ay exact caleuletron with drag might be consrolered wm 
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Range correction 


5 Cy ine wy dt AL 
(33) A pate zp f/t> 


and for the sideward deviations caused by a wind sae? w, 
perpendicular to the missile velocity, v, 


(3h) ee ee sel]? ‘ar We dt 


Since the trajectory is known, the integrals (33) and (34) can 
easily be evaluated. In many cases it may be sufficient to use 
an average speed, v_, between burnout and the peak and to use for 


the variable horizontal speed, v,, the horizontal speed at burnout, 
Vp}> Then the equations are simplified as (35) and (36); by ty 
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Fig. 4. COORDINATES ON THE SPHERE THROUGH BURNOUT 
POINT B AND TARGET POINT T 
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